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Abstract. The aim of this paper is to consider quasi- linear systems 
which are not in the form of evolution equations. We propose new 
condition of Richness or Semi-Hamiltonicity for such a system and prove 
that the blow up analysis along characteristic curves can be performed 
for it in an analogous manner. This opens a possibility to use this ansatz 
also for geometric problems. We apply the results to the problem of 
Polynomial integral for geodesic flows on the 2-torus. 



1. Motivation and the result 

Consider a quasi-linear system for vector function u(x,y) = (ui,...,u n ) 
which has the following form 

A(u)u x + B(u)u y = (1) 

It may happen in practice that one of the matrices A(u) and B(u) can 
degenerate somewhere (and even both of them can degenerate somewhere). 

Our main assumption about these matrices is that the homogeneous poly- 
nomial P in a, (3 is not a zero polynomial at any point (x, y) : 

P = det(aB - /3A), deg(P) = n. (P) 

This assumption is obviously satisfied if one of the matrices A(u) and B(u) 
is non-degenerate, however we shall assume everywhere the weaker version- 
(P). We shall see in the example that (P) is in fact the correct assumption. 

Moreover, we shall assume in the following that the system is strictly 
hyperbolic that is the polynomial P has n distinct roots [/% : a,]. We define 
unite characteristic vector fields on the plane R 2 (x,y) by 



Vi = cos 4>id x + sin 4>id y , 

where the angles (pi, we shall call them characteristic angles, are such that 
[sin^j : cos<^,] = [ft : oij]. 

Notice that if one of the matrices, say A is non-degenerate then the system 
is equivalent to one in the evolution form. The notion of the system to be 
Rich or Semi-Hamiltonian (see |1U|.|11] and |12j.[5], we shall call them Rich 
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for the sake of brevity) for evolution system says that it can be written in 
Riemann invariants (diagonal form) 

(ri) x + Aj(n, ...,r n )(ri) y = 0,i = 1, ...,n, 

and moreover the eigenvalues Aj = (3i/cti of A~~ 1 B satisfy the following 
identities: 

/ d n Xj \ = ^ / d rk Xj \ _ ^ 

rfc \Xi - XjJ Ti \x k - XjJ 

This condition allows one to perform blow-up analysis along characteristics 
as it is shown in |10| and used for example in pjj. It was proved by B. 
Sevennec [IT] and later understood in differential-geometric terms [7] that 
strictly Hyperbolic system in evolution form which is written in Riemann 
invariants is Rich, if and only if there are local coordinates in which the 
system takes the form of conservation laws. 

The unsatisfactory thing, however, with the condition (R) is the fact that 
for the system (1) characteristic curves can pass from the chart where A 
is non-degenerate to the chart where B is non-degenerate or even reach 
those points where both matrices degenerate. This does-not allow to use 
the analysis of corresponding Riccati equations for all times. 

We propose the following generalization of the Richness condition whose 
naturality we shall justify below: 

Definition 1.1. We call system (1) Rich if it can be written in the diagonal 
form 

L Vi n = cos <j>i(ri) x + sn\(j)i(r i )y = 0, i = l,...,n (2) 

for a regular change of variables (u%, ...,u n ) — > (r\, r n ) and the following 
conditions on the characteristic angles <f>i(ri, ...,r n ) holds true 

° Tk Vtan(^ - 4>j)J n Vtan(«^ - fa) J { ] 

It is important fact that this definition is invariant with respect to rota- 
tions of the plane. We shall continue to call r{ in (2) by Riemann invariants. 
Our first result is the following 

Theorem 1.2. If system (1) satisfying (P) is Rich according Definition 1.1, 
so that the conditions (2),(<&) hold true, then the derivatives of i-th Riemann 
invariant Wi = L ±n in the orthogonal direction to characteristics satisfy 

i 

the following Riccati equation: 

L^(exp {-Gi)wi) + exp (Gi)d ri ((pi)(exp(-Gi)wi) 2 = 0, 
for any i = 1, ...,n, where Gj is a function of Riemann invariants satisfying 

r 



Or j G 7 



3 tan(<£i - fa)' 

Here v 1 - stands for the vector field rotated from v by 90° counterclockwise. 

We shall see in lemma below that the conditions (R) and ($) are almost 
equivalent. This lemma enables us to prove the following theorem which is 
a generalization to our case of the result of 
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Theorem 1.3. Given any strictly Hyperbolic diagonal system 

cos <pi(ri) x + sin 0j (r^ = 0, i = 1, ...,n. 

Then the condition (<&) is satisfied if and only if the system can be written 
in the form of n conservation laws 

{9i)x + {hi) y = 0,i = l,...,n. 

Remark 1. Saying that system (1) can be written in a certain form, one 
means by definition that this form can be achieved by a change of variables 
and multiplication by an invertible matrix: Let u = <3?(u?) be a regular 
change of variables and C(w) is an invertible matrix then one can write 
system (1) in the form 

C(w)A($(w))D<S>(w)w x + C(w)B($(w))D<S>(w)w y = 0, 

where D is the differential. 
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2. Derivation along characteristics. Proof of Theorem 1.2. 
Differentiate the j-th equation of (2) with respect to the field vj~. We have 
= L v ± L Vj rj = L Vj L v ± rj - L [v . tV ± } rj . (3) 
Compute now the derivative along the commutator: 

L [vj,vf] r 3 = L Vj L vf r i ~ L vf L Vj r j = 

= L Vj (- siD.(f>j(rj) x + cos<j)j(rj)y) - L v ± (cos <pj(rj) x + sin <pj(rj) y ) = 

= {fj) x {- COS 2 <f>j((f>j) x -COS(f>j8Ul(/>j(<f>j)y) + 

+ (rj) y (- sin. <f>j cos <f>j(<f>j) x - sin 2 <f>j((j>j)y)+ 

+ { r j)x{- sin 2 <pj{<j)j) x + sin^ cos 4>j{4>j) y ) + 

+ (r j) y (sin <j) j cos 4>j(<j)j) x - cos 2 (f>j{(f>j)y) = 

= -(rj)x(<Pj)x ~ (rj)y(<Pj) y (4) 
Notice that the derivatives (rj) x , (rj) y can be expressed by the following 
two identities: 

cos cf>j(rj) x + sin<j)j(rj)y = 0, 
- sin^-(rj)x + cos^-fa),, = L v ±rj. 

Therefore 

{tj) x = -am<i>jL v ±r j , 

(rj) x = cos <PjL v ± rj . (5) 
Substituting back to (4) we get 

L [v 3 ,v^\ r i = {L v *-r j )(wa.<l>j{<l>j) x - CQBfatyrfy) = -(L v ±rj)(L v ±<f>j) (6) 
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By the chain rule for L v ±<f>j the last equation can be rewritten as follows 

n 

L [« i) «+] r J = - L vf r j^^ d ri^ L vf r i) = 
3 3 i=l J 

-{L^Tjfidr^j) - {L vf r 3 ) Y,{dr^ 3 ){L vf n). (7) 

Let me express now the derivative 

L v ±n = - sm(/) j (r i ) x + cos(f>j(ri) y , (8) 

via L Vj ri as follows. Write 

cos^j(rj)a; + sin 4>i(ri) y = 0, 

cos0j(rj)a; + sin 4>j(ri)y = L Vj r,i. (9) 
From these two identities we have 

smOi - 0j) 3 s\n(4>i -4>j) 3 

Substitute expressions (10) into (8) to get: 

Ly - Ti 



Denote by 



Wi := L v ±n. 



Plug this together with (11) into equation (7) and then to (3): 

L Vj (wj) + (dr^yiwj) 2 - Wj Y^( dr ih \ an u- - 6 ) Lv i ri = °' ^ 



By Richness ($>) we have that for all j = 1, n there exist functions 

0j (r,,...,r„): *,e i = i -^L J , j ^. (13) 
By (13) we can rewrite (12) as 



£ty(u>j) + {d ri <t>j)(wj) 2 - Wj Y^{driGj)L Vj ri = 0, 

which is the same as 

L Vj ( Wj ) + (dr^ywj) 2 - wjL^Gj = 0. (14) 

Multiplying (14) by exp(-Gj) we get the Riccati equation of the first theo- 
rem for j instead of i. This completes the proof of theorem 1.2. 
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3. Conservation laws. Proof of theorem 1.3. 

We shall need the following key observation. 

Lemma 3.1. Given two sets of functions Aj(ri, ...,r n ); <f>i(r±, ■■■r n ), i = 
1, n) such that 

Aj Xj, (pi it/2 (modir), Xi = tan^. 

Then the conditions (R) and (*f>) are equivalent. 

The proof which I know is computational. It would be interesting to find 
more conceptual proof. 

Proof. Let us prove first that (R) implies (<£). 
Denote by 

3 n Xj 



J Xi — Xj 
Then satisfy the following identities (|10j): 

dri®kj = drf.O'ij = akiaij + aik&kj &kjQ'ij- (15) 

In order to prove them differentiate with respect to the identity 

d n Xj = aij(Xi — Xj) then interchange the order of i, k, subtract one from 

the other and divide by Aj — A&. 
Denote by 



d ri <t>j _ d n Xj 1 + XjXj _ 1 + XjXj 
bij - tan(0 4 - 4>j) ~ Xi - X, 1 + A] " ^ 1 + X] ' 



To prove ($) we have to verify that the difference 

d = d Tk bij — d n bkj 
vanishes. Let us compute d explicitly: 



7 /rj N l+AjAj 1 + XkXj 

d = (d^aij) 1 + A 2 - (^-iOfcj) 1 + A 2 



.7 

, „ / 1 + AiAjA . / 1 + X k Xj\ 

+« A l^TTAf ) " flfcAl I^TTAf J ■ 

By the identities (15) and the condition (R) we have 

, / . \ Aj(Aj — Afe) 

" = [akiaij + aikOkj — akjaij) — 2 h 

l + \j 

a rfc (A 4 A i )(l + A|) - (1 + AjAj)2Aj9 rfc (Aj) 

+ CLij 



(1 + A?) 2 



a r ,(A fc A i )(l + A| ) - (1 + X k Xj)2Xjd ri (Xj) 



n ! \5\2 • ( i7 ) 

V 

Substitute now into the nominators of (17) the following expressions for 



the derivatives of Xj from the definition of a^: 



d Ti Xj — dijiXi — Xj). 
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Then one has 



, Aj(Aj ~ Afc) 



d — (flki&ij ~\~ O'ikO'kj 0, k jdij) i i \2 

i + A i 

(a k j(X k - Xj)Xj + Xja ki (X k - Xi) (aijjXj - Xj)X k + XjOjk(Xi - X k ) 
+<Hj 1 + A 2 a kj 1 + A 2 

(1 + XiXj)Xja kj (X k - Xj) (1 + A/,A / ;A J (/, ; iA, - Xj) 

(TTaJP + 2afeJ oTa|? ' ( } 

Notice that the identity (18) is a quadratic expression in a^s. Collecting 
the coefficients of aija k j, ai k a k j, a k idij one comes to d = 0. This proves 
lemma in one direction. 

Proof of the converse statement is very much analogous but with even 
harder computations. I shall reproduce them sketchy. So we assume the 
identities (<£) are satisfied. First one can obtain the identity analogous to 
(15) for the derivatives d Tk bij in the following way. Write 

d n (pj = hjt&n^i - <j>j) = bij ^ , d n Xj = fr / + 1 ^* 3 \ (19). 

Differentiating the first equality of (19) with respect to r k , using the identi- 
ties (19) again and taking into account (16) one has 



d rk d ri <f>j = d rk {bij)- 



(Aj — Xj) 



+ AjA, 



+ 



•j 



(Aj — Xj) 2 \ ( X k — Xi X k — X 



(1 + XiXj) 2 J V 1 + AfeAi KJ 1 + AfeAj-, 
Interchanging in this identity the order of indexes i and k and using d rk bij 
d ri b k j one has the identity: 

. , (A,-A fc )(l + A 2 ) 
^^ j (l + A l A,)(l + A fe A,) " 

-h (a 4. ( A fc ~ A i) 2 A / . A » ~ A fc . A * - A i 



(1 + WA 1 + AiAfc J 1 + AjAj 



(1 + AiAj) 2 / V 1 + AiA fe 1 + AfeAj 
In order to verify (R) one computes 

A,(A fc -A,)(l + A 2 ) 



^ r k ^kj — &ij 



(1 + XiXj)(l + X k Xj 



+i 'A(i^)- 6 'A(iT^-)- < 21 » 

The last step is to plug into (21) the expression (20) and also to differentiate 
the last two brackets of (21) using the expression for the derivatives (19). 
Then one finally gets a quadratic expression in b%jS. Collecting similar terms 
one verifies that the right hand side of (21) vanishes. Therefore (R) holds 
true. This proves the lemma. □ 

It is easy now to prove Theorem 1.3. 
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Proof. Notice first of all that the statement of the second theorem is local. 
Given a system which is strictly Hyperbolic and is written in the diagonal 
form 

cos 4>i(ri) x + sin (f>i(ri) y = 0, i = 1, ...,n, 

Let us give a proof first in one direction, namely assume that the system 
can be written in the form of conservation laws 

(9i)x + (hi) y = 0, i = 1, ...,n. 

Let me explain that then it must satisfy condition (<&). If among fa there is 
one with cos fa = then one can apply a small rotation of the plane R 2 (x, y) 
and to get a new system which has all angles different from ±7r/2. Notice 
that the rotated system remains in the form of conservation laws and in 
addition the differential Dg becomes a non singular matrix, since otherwise 
a = 1, (3 = would be the root of (P) but this is impossible by fa / ±7r/2. 
Denote 

Aj := tan fa. 

Use now Sevennec' theorem saying that the diagonal system 

(n) x + Xi(n) y = o 

which can be written in the form of conservation laws 

(9i)x + (hi) y = 0,i = 1, ...,n 

with the non-singular Jacobi matrix (d rj gi) must satisfy (R). But by lemma 
in this case (R) and ($) are equivalent. So we get condition ($) for rotated 
system. But this condition is obviously rotationally invariant. Thus it holds 
also for the original system. 

The proof in the opposite direction is analogous. First rotate the plane 
exactly as above. Condition (<£)remains valid since it is rotationally invari- 
ant. Then by the lemma (R) is valid as well and then by Sevennec' theorem 
the rotated system can written in the form of conservation laws. But then 
obviously the original one as well. This completes the proof. 

□ 

4. Geometric example. 

In this section we give a geometric example where the results of the pre- 
vious sections are important. 

Let p be a Riemannian metric on the 2-torus T 2 = R 2 /T, p* denotes the 
geodesic flow. Assume that p is written in conformal way: 

ds 2 = A(x,y)(dx 2 + dy 2 ). 

Let F : T*T 2 be a function on the cotangent bundle which is homogeneous 
polynomial of degree n with respect to the fibre: 

n 
k=0 

We are looking for such an F which is an integral of motion for the geodesic 
flow p t , i.e. Fop 1 = F. We shall also assume that this F is irreducible, i.e of 
minimal possible degree. Let us mention that this problem is classical there 
are very well studied examples of the geodesic flows on the 2-torus which 
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have integrals F of degree one and two. We refer to books [I] and [9] for 
the history and discussion of this classical question with references therein. 
In our recent papers with A.E.Mironov we used the so called semi-geodesic 
coordinates. In these coordinates one arrives to a remarkable Rich quasi- 
linear system of equations in evolution form on the coefficients of the integral 
F (|2]> ED- However it is very natural to be able to work in conformal 
coordinates as well. In this case the quasi-linear system on the coefficients 
has no evolution form any more but looks like: 

A(U)U x + B(U)U y = 0. 

Let me write down explicitly the matrices for the case n = 3 (this case is 
already very interesting and not trivial see for example [6]). 

1 3a 

A(U) =| 1 36 | , B{U) = | 1 -3a | , U = | v | . (22) 
A u 

Here a, 6, u, v are related to the coefficients of the integral a, by the following 

U V U V 

a = a + — , ai = 36 + — , a 2 = -3a + — , a 3 = -6 + — . 

It was noticed in [S] that a, 6 are in fact constants. Computing polynomial 
P = det(aB — (3 A) one has: 

P = a 3 (v + 36A) + a 2 f3(-u - 9aA) + a/3 2 (v - 96A) + (3 3 (-u + 3aA). 

Let remark that it may happen at some points that both matrices A, B are 
degenerate, however polynomial P for any point can not vanish identically. 
This is because otherwise both constants a, 6 vanish, but then one checks 
that in such a case the integral F is a product of the Hamiltonian with an 
integral of degree one in momenta, therefore reducible. 

Notice that quasi-linear system (22) is written in the form of conservation 
laws 

{gi)x + (hi) y = 0, 

<?i = u + 3aA, 02 = v + 36A, 53 = uA, 

h\ = —v + 36A, h% = u — 3aA, /13 = vA. 

Moreover by a very general argument in the Hyperbolic region this system 
can be written in the diagonal form (2). Indeed introduce angular coordinate 
<f) on the fibres of the energy level 

{t(Pi +P2) = 1} : Pi = \/Acos<?!>,p 2 = \/Asin0, 

then one can verify that the condition on a function F to be an integral of 
the flow reads 



2A T 2A 

At the points where Fa, vanishes this equation takes particularly nice form: 

F x cos 4> + F y sin = 0. 

Therefore critical values of F on the fibre are Riemann invariants. One can 
check also that the polynomial P is proportional in fact to the derivative 
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of F in the direction of the fibre. Moreover one can check, as we did in 
[2], that in the Hyperbolic region Riemann invariants form a regular change 
of variables. As a consequence of Theorem 1.3 one concludes that in the 
Hyperbolic region the system of this example is Rich in our generalized 
sense. And therefore Theorem 1.2 tells us that the Riccati equation along 
characteristics applies. 

5. Questions 
Several questions are very natural: 

1. It would be interesting to find more conceptual proof of the lemma in 
the framework differential-geometric approach by Dubrovin-Novikov [5]. 

2. How does generalized Hodograph method by Tsarev [12] work in our 

case? 

3. How to analyze the behavior of the Riccati equation for the example 
of previous section? It seems that genuine non-linearity condition can not 
be expected for all eigenvalues. 
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